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The  constant  energy  dissipative  particle  dynamics  (DPD-E)  method  is  implemented  into  the 
Large-Scale  Atomic/Molecular  Massively  Parallel  Simulator  (LAMMPS)  simulation  software  to 
efficiently  model  systems  under  isoenergetic  conditions  using  the  current  LAMMPS  velocity- 
Verlet  (VV)  integration  scheme  and  the  Shardlow-splitting  algorithm  (SSA).  The  relevant  source 
code  files  that  are  current  with  the  7  Lebruary  2014  release  of  LAMMPS  are  provided,  along 
with  the  user  manual  documentation.  The  contents  of  this  material  are  described  in  detail  in  U.S. 
Army  Research  Laboratory  (ARL)  technical  report  ARL-TR-6863.1 

The  DPD-E  method  using  the  VV  and  VV-SSA  integration  schemes  has  been  implemented 
within  LAMMPS  as  user  packages  under  the  directory  src/USER-DPD  in  the  accompanying 
compact  disc  (CD). 

The  LAMMPS  user  package  USER-DPD  can  be  added  to  the  LAMMPS  source  files  by  issuing 
the  command 

make  yes-USER-DPD 

within  the  LAMMPS  src/  directory.  A  copy  of  the  modified  user  manual  containing  the  new 
features  added  to  LAMMPS  can  be  found  in  the  appendix,  and  the  source  txt,  html,  tex,  and  jpg 
files  can  be  found  within  the  LAMMPS  doc/  and  doc/Eqs  directory  in  the  accompanying  CD. 

The  full  LAMMPS  user  manual  can  be  found  at  http://lammps.sandia.gov/doc/Manual.html. 


1  Larentzos,  J.  P.;  Brennan,  J.  K.;  Moore.  J.  D.;  Mattson,  W.  D.  LAMMPS  Implementation  of  Constant  Energy  Dissipative 
Particle  Dynamics  (DPE-E);  ARL-TR-6863;  U.S.  Army  Research  Laboratory:  Aberdeen  Proving  Ground,  March  2014. 
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Appendix.  Large-Scale  Atomic/Molecular  Massively  Parallel  Simulator 
(LAMMPS)  User  Manual  for  Commands  Related  to  the  Constant  Energy 
Dissipative  Particle  Dynamics  (DPD-E)  Implementation  With  the  Velocity  - 
Verlet  (VV)  and  VV-Shardlow-Splitting  Algorithm  (SSA)  Integration 

Schemes 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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atom  style  command 

Syntax: 

atom_style  style  args 

•  style  =  angle  or  atomic  or  body  or  bond  or  charge  or  dpd  or  dipole  or  electron  or  ellipsoid  or  full  or 
line  or  meso  or  molecular  or  peri  or  sphere  or  tri.  or  hybrid 

args  =  none  for  any  style  except  body  and  hybrid 
body  args  =  bstyle  bstyle-args 
bstyle  =  style  of  body  particles 

bstyle-args  =  additional  arguments  specific  to  the  bstyle 
see  the  body  doc  page  for  details 

hybrid  args  =  list  of  one  or  more  sub-styles,  each  with  their  args 


Examples: 


atom_style 

atom_style 

atom_style 

atom_style 

atom_style 

atom_style 


atomic 

bond 

full 

body  nparticle  2  10 

hybrid  charge  bond 

hybrid  charge  body  nparticle  2  5 


Description: 


Define  what  style  of  atoms  to  use  in  a  simulation.  This  determines  what  attributes  are  associated  with  the 
atoms.  This  command  must  be  used  before  a  simulation  is  setup  via  a  read  data,  read  restart,  or  create  box 
command. 


Once  a  style  is  assigned,  it  cannot  be  changed,  so  use  a  style  general  enough  to  encompass  all  attributes.  E.g. 
with  style  bond ,  angular  terms  cannot  be  used  or  added  later  to  the  model.  It  is  OK  to  use  a  style  more  general 
than  needed,  though  it  may  be  slightly  inefficient. 

The  choice  of  style  affects  what  quantities  are  stored  by  each  atom,  what  quantities  are  communicated 
between  processors  to  enable  forces  to  be  computed,  and  what  quantities  are  listed  in  the  data  file  read  by  the 
read  data  command. 

These  are  the  additional  attributes  of  each  style  and  the  typical  kinds  of  physical  systems  they  are  used  to 
model.  All  styles  store  coordinates,  velocities,  atom  IDs  and  types.  See  the  read  data,  create  atoms,  and  set. 
commands  for  info  on  how  to  set  these  various  quantities. 


angle 

bonds  and  angles 

bead-spring  polymers  with  stiffness 

atomic 

only  the  default  values 

coarse-grain  liquids,  solids,  metals 

body 

mass,  inertia  moments,  quaternion,  angular  momentum 

arbitrary  bodies 

bond 

bonds 

bead-spring  polymers 

charge 

charge 

atomic  system  with  charges 

dpd 

internal  temperature  and  energies 

DPD  particles 

atom_style  command 
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dipole 

charge  and  dipole  moment 

system  with  dipolar  particles 

electron 

charge  and  spin  and  eradius 

electronic  force  field 

ellipsoid. 

shape,  quaternion,  angular  momentum 

aspherical  particles 

mmm 

molecular  +  charge 

bio-molecules 

line 

end  points,  angular  velocity 

rigid  bodies 

meso 

rho,  e,  cv 

SPH  particles 

molecular 

bonds,  angles,  dihedrals,  impropers 

uncharged  molecules 

peri 

mass,  volume 

mesocopic  Peridynamic  models 

sphere 

diameter,  mass,  angular  velocity 

granular  models 

tri 

comer  points,  angular  momentum 

rigid  bodies 

wavepacket 

charge,  spin,  eradius,  etag,  cs_re,  cs_im 

AWPMD 

All  of  the  styles  define  point  particles,  except  the  sphere,  ellipsoid,  electron,  peri,  wavepacket,  line,  tri,  and 
body  styles,  which  define  finite-size  particles.  See  Section  howto  14  for  an  overview  of  using  finite-size 
particle  models  with  LAMMPS. 


All  of  the  styles  assign  mass  to  particles  on  a  per-type  basis,  using  the  mass  command,  except  for  the 
finite-size  particle  styles.  They  assign  mass  to  individual  particles  on  a  per-particle  basis. 

For  the  sphere  style,  the  particles  are  spheres  and  each  stores  a  per-particle  diameter  and  mass.  If  the  diameter 
>  0,0,  the  particle  is  a  finite-size  sphere.  If  the  diameter  =  0.0,  it  is  a  point  particle. 

For  the  ellipsoid  style,  the  particles  are  ellipsoids  and  each  stores  a  flag  which  indicates  whether  it  is  a 
finite-size  ellipsoid  or  a  point  particle.  If  it  is  an  ellipsoid,  it  also  stores  a  shape  vector  with  the  3  diamters  of 
the  ellipsoid  and  a  quaternion  4- vector  with  its  orientation. 

For  the  electron  style,  the  particles  representing  electrons  are  3d  Gaussians  with  a  specified  position  and 
bandwidth  or  uncertainty  in  position,  which  is  represented  by  the  eradius  =  electron  size. 

For  the  peri  style,  the  particles  are  spherical  and  each  stores  a  per-particle  mass  and  volume. 

The  meso  style  is  for  smoothed  particle  hydrodynamics  (SPH)  particles  which  store  a  density  (rho),  energy 
(e),  and  heat  capacity  (cv). 

The  dpd  style  is  for  dissipative  particle  dynamics  (DPD)  particles  which  store  the  particle  internal  temperature 
(dpdTheta),  internal  conductive  energy  (uCond),  internal  mechanical  energy  (uMech),  as  well  as  other  particle 
properties  such  as  heat  capacity  (cv)  and  density  (rho). 

The  wavepacket  style  is  similar  to  electron,  but  the  electrons  may  consist  of  several  Gaussian  wave  packets, 
summed  up  with  coefficients  cs=  (cs_re,cs_im).  Each  of  the  wave  packets  is  treated  as  a  separate  particle  in 
LAMMPS,  wave  packets  belonging  to  the  same  electron  must  have  identical  etag  values. 

For  the  line  style,  the  particles  are  idealized  line  segments  and  each  stores  a  per-particle  mass  and  length  and 
orientation  (i.e.  the  end  points  of  the  line  segment). 

For  the  tri  style,  the  particles  are  planar  triangles  and  each  stores  a  per-particle  mass  and  size  and  orientation 
(i.e.  the  comer  points  of  the  triangle). 
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For  the  body  style,  the  particles  are  arbitrary  bodies  with  internal  attributes  defined  by  the  "style”  of  the 
bodies,  which  is  specified  by  the  bstyle  argument.  Body  particles  can  represent  complex  entities,  such  as 
surface  meshes  of  discrete  points,  collections  of  sub-particles,  deformable  objects,  etc. 

The  body  doc  page  descibes  the  body  styles  LAMMPS  currently  supports,  and  provides  more  details  as  to  the 
kind  of  body  particles  they  represent.  For  all  styles,  each  body  particle  stores  moments  of  inertia  and  a 
quaternion  4-vector,  so  that  its  orientation  and  position  can  be  time  integrated  due  to  forces  and  torques. 

Note  that  there  may  be  additional  arguments  required  along  with  the  bstyle  specification,  in  the  atom  style 
body  command.  These  arguments  are  described  in  the  body  doc  page. 


Typically,  simulations  require  only  a  single  (non-hybrid)  atom  style.  If  some  atoms  in  the  simulation  do  not 
have  all  the  properties  defined  by  a  particular  style,  use  the  simplest  style  that  defines  all  the  needed  properties 
by  any  atom.  For  example,  if  some  atoms  in  a  simulation  are  charged,  but  others  are  not,  use  the  charge  style. 
If  some  atoms  have  bonds,  but  others  do  not,  use  the  bond  style. 

The  only  scenario  where  the  hybrid  style  is  needed  is  if  there  is  no  single  style  which  defmes  all  needed 
properties  of  all  atoms.  For  example,  if  you  want  dipolar  particles  which  will  rotate  due  to  torque,  you  would 
need  to  use  "atom  style  hybrid  sphere  dipole".  When  a  hybrid  style  is  used,  atoms  store  and  communicate  the 
union  of  all  quantities  implied  by  the  individual  styles. 

LAMMPS  can  be  extended  with  new  atom  styles  as  well  as  new  body  styles;  see  Ibis  section. 

Restrictions: 

This  command  cannot  be  used  after  the  simulation  box  is  defined  by  a  read  data  or  create  box  command. 

The  angle ,  bond,  full,  and  molecular  styles  are  part  of  die  MOLECULAR  package.  The  line  and  tri  styles  are 
part  of  the  ASPHERE  pacakge.  The  body  style  is  part  of  the  BODY  package.  The  dipole  style  is  part  of  the 
DIPOLE  package.  The  peri  style  is  part  of  the  PERI  package  for  Peridynamics.  The  electron  style  is  part  of 
the  USER-EFF  package  for  electronic  force  fields.  The  dpd  style  is  part  of  the  USER-DPDE  and 
USER-DPDE-SHARDLOW  packages  for  dissipative  particle  dynamics  (DPD).  The  meso  style  is  part  of  the 
USER-SPH  package  for  smoothed  particle  hydrodyanmics  (SPH).  See  this  PDF  guide  to  using  SPH  in 
LAMMPS.  The  wavepacket  style  is  part  of  the  USER-AWPMD  package  for  the  anti  symmetrized  wave  packet 
MD  melhod.  They  are  only  enabled  if  LAMMPS  was  built  with  that  package.  See  the  Making  LAMMPS 
section  for  more  info. 

Related  commands: 

read  data,  pair  style 

Default: 

atom  style  atomic 


atom_style  command 


3 


6 


LAMMPS  WWW  Site  -  LAMMPS  Documentation  -  LAMMPS  Commands 


compute  dpd  command 

Syntax: 

compute  ID  group- ID  dpd 

•  ID,  group -ID  are  documented  in  compute  command 

•  dpd  =  style  name  of  this  compute  command 

Examples: 

compute  1  ell  dpd 

Description: 


Define  a  computation  that  accumulates  the  total  internal  conductive  energy  (U_cond) ,  the  total  internal 
mechanical  energy  (U_mech),  the  total  internal  energy  (U)  and  the  average  internal  temperature  (Theta)  of  the 
entire  system  of  particles.  See  the  compute  dp  d/atom  command  if  you  want  per-particle  internal  energies  and 
internal  temperatures. 


The  system  internal  properties  are  computed  according  to  the  following  relations: 

N 


jjcond  =  ^2  u': 
i= 1 
N 

jjmech  _  ^  ^  yjji 


cond 


rnech 


i=l 


N 


u  =  Y,Kond  +  u?ech) 


2=1 


i  A  i 


where  N  is  the  number  of  particles  in  the  system 


Output  info: 

This  compute  calculates  a  global  vector  of  length  4  (U_cond,  U_mech,  U,  Theta),  which  can  be  accessed  by 
indices  1-4.  See  this  section  for  an  overview  of  LAMMPS  output  options. 

The  vector  values  will  be  in  energy  and  temperature  units. 

Restrictions: 

The  compute  dpd  is  only  available  if  LAMMPS  is  built  with  the  appropriate  USER-DPDE  or 
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USER-DPDE-SHARDLOW  package. 
Related  commands: 
compute  dtxl/atoin.  thermo  style 
Default:  none 


(Larentzos)  J.P.  Larentzos,  J.K.  Brennan,  J.D.  Moore,  and  W.D.  Mattson,  "LAMMPS  Implementation  of 
Constant  Energy  Dissipative  Particle  Dynamics  (DPD-E)”,  ARL-TR-XXXX,  U.S,  Army  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD, 
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compute  dpd/atom  command 

Syntax: 

compute  ID  group- ID  dpd/atom 


•  ID,  group-ID  are  documented  in  compute  command 

•  dpd/atom  =  style  name  of  this  compute  command 

Examples: 

compute  1  all  dpd/atom 
Description: 

Define  a  computation  that  accesses  the  per-particle  internal  conductive  energy  (u  cond),  internal  mechanical 
energy  (u_mech)  and  internal  temperatures  (  theta)  for  each  particle  in  a  group.  See  the  compute  dpd  command 
if  you  want  the  total  internal  conductive  energy,  the  total  internal  mechanical  energy,  and  average  internal 
temperature  of  the  entire  system  of  dpd  particles. 

Output  info: 

This  compute  calculates  a  per-particle  array  with  3  columns,  which  can  be  accessed  by  indices  1-3  by  any 
command  that  uses  per-particle  values  from  a  compute  as  input.  See  Seel  ion  howto  15  for  an  overview  of 
LAMMPS  output  options. 

The  per-particle  array  values  will  be  in  energy  and  temperature  units  as  discussed  above. 

Restrictions: 

The  compute  dpd  is  only  available  if  LAMMPS  is  built  with  the  appropriate  USER-DPDE  or 
USER-DPDE-SHARDLOW  package. 

Related  commands: 

compute  dpd 

Default:  none 


(Larentzos)  J.P.  Larentzos,  J.K.  Brennan,  J.D.  Moore,  and  W.D.  Mattson,  "LAMMPS  Implementation  of 
Constant  Energy  Dissipative  Particle  Dynamics  (DPD-E)”,  ARL-TR-XXXX,  U.S.  Army  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD. 
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fix  dpde  command 

Syntax: 


fix  ID  group-ID  dpde 


ID,  group-ID  are  documented  in  fix  command  dpde  =  style  name  of  this  fix  command 
Examples: 


fix  1  all  dpde 


Description: 

Perform  constant  energy  dissipative  particle  dynamics  (DPD-E)  integration  to  update  position,  velocity  and 
internal  energy  for  particles  in  the  group  at  each  timestep. 

For  fix  dpde ,  the  particle  internal  temperature  is  related  to  the  particle  internal  energy  through  a  mesoparticle 
equation  of  state.  An  additional  fix  must  be  specified  that  defines  the  equation  of  state  for  each  particle. 


Restrictions: 

The  fix  dpde  is  only  available  if  LAMMPS  is  built  with  the  appropriate  USER-DPDE  package. 

The  fix  dpde  requires  the  dpd  atom  style  to  be  used  in  order  to  properly  account  for  the  particle  internal 
energies  and  temperature. 

The  fix  dpde  must  be  used  with  an  additional  fix  that  specifies  the  mesoparticle  equation  of  state  for  each 
particle. 

Related  commands: 

fix  live  fix  eos/cv 

Default:  none 


(Larentzos)  J.P.  Larentzos,  J.K.  Brennan,  J.D.  Moore,  and  W.D.  Mattson,  "LAMMPS  Implementation  of 
Constant  Energy  Dissipative  Particle  Dynamics  (DPD-E)",  ARL-TR-XXXX,  U.S.  Army  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD. 
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fix  dpde/shardlow  command 

Syntax: 


fix  ID  group-ID  dpde/shardlow 


ID,  group-ID  are  documented  in  fix  command  dpde/shardlow  =  style  name  of  this  fix  command 
Examples: 

fix  1  all  dpde/ shardlow 

Description: 

Perform  constant  energy  dissipative  particle  dynamics  (DPD-E)  integration  using  the  Shardlow  splitting 
algorithm  (SSA)  to  update  position,  velocity  and  internal  energy  for  particles  in  the  group  at  each  timestep. 
The  SSA  splits  the  integration  into  a  stochastic  and  deterministic  integration  step.  The  stochastic  integration 
of  the  dissipative  and  random  forces  is  performed  prior  to  the  deterministic  integration  of  the  conservative 
force.  Further  details  regarding  the  method  are  provided  in  (T.isal)  and  (T.arentzos). 

For  fix  dpde/shardlow,  the  particle  internal  temperature  is  related  to  the  particle  internal  energy  through  a 
mesoparticle  equation  of  state.  An  additional  fix  must  be  specified  that  defines  the  equation  of  state  for  each 
particle. 


Restrictions: 

The  fix  dpde/shardlow  is  only  available  if  LAMMPS  is  built  with  the  appropriate 
USER-DPDE-SHARDLOW  package. 

The  fix  dpde/shardlow  must  be  used  with  the  dpd/conservative  pair  style  command  to  properly  initialize  pair 
coefficients  for  sigma  and  kappa. 

The  fix  dpde/shardlow  requires  the  dpd  atom  style  to  be  used  in  order  to  properly  account  for  the  particle 
internal  energies  and  temperature. 

The  fix  dpde/shardlow  must  be  used  with  an  additional  fix  that  specifies  the  mesoparticle  equation  of  state  for 
each  particle. 

Related  commands: 

fix  dnde  fix  eos/cv 

Default:  none 


(Lisal)  M.  Lisal,  J.K.  Brennan,  J.  Bonet  Avalos,  Dissipative  particle  dynamics  as  isothermal,  isobaric, 
isoenergetic,  and  isoenthalpic  conditions  using  Shardlow-like  splitting  algorithms.”,  I.  Chem.  Phys.,  135, 
204105(2011). 
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(Larentzos)  J.P.  Larentzos,  J.K.  Brennan,  J.D.  Moore,  M.  Lisal  and  W.D.  Mattson,  "Parallel  Implementation 
of  Isothermal  and  Isoenergetic  Dissipative  Particle  Dynamics  Using  Shardlow-Like  Splitting  Algorithms”, 
Submitted  to  CompuL  Phys,  Commun.,  (2013). 


(Larentzos)  J.P.  Larentzos,  J.K.  Brennan,  J.D.  Moore,  and  W.D.  Mattson,  "LAMMPS  Implementation  of 
Constant  Energy  Dissipative  Particle  Dynamics  (DPD-E)'1,  ARL-TR-XXXX,  U.S.  Army  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD. 
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fix  eos/cv  command 

Syntax: 

fix  ID  group-ID  eos/cv  cv 

•  ID,  group-ID  are  documented  in  fix  command 

•  eos/cv  =  style  name  of  this  fix  command 

•  cv  =  constant-volume  heat  capacity 

Examples: 

fix  1  all  eos/cv  0.01 

Description: 

Fix  eos/cv  applies  a  mesoparticle  equation  of  state  to  relate  the  particle  internal  energy  (u_i)  to  the  particle 
internal  temperature  ( theta  i).  The  eos/cv  mesoparticle  equation  of  state  requires  the  constant-volume  heat 
capacity,  and  is  defined  as  follows: 


,  ,.cond  /-<  a 
Ui  =Ut  +  U{  CV0i 

where  Cv  is  the  constant- volume  heat  capacity,  u_cond  is  the  internal  conductive  energy,  and  u_mech  is  the 
internal  mechanical  energy.  Note  that  alternative  definitions  of  the  mesoparticle  equation  of  state  are  possible, 
but  not  currently  implemented. 


Restrictions: 

The  fix  eos/cv  is  only  available  if  LAMMPS  is  built  with  the  appropriate  USER-DPDE  or 
USER-DPDE-SHARDLOW  packages. 

The  fix  eos/cv  must  be  used  with  the  atom  style  dpd. 

Related  commands: 
fix  dpde 
Default:  none 


(Larentzos)  J.P.  Larentzos,  J.K.  Brennan,  J.D.  Moore,  and  W.D.  Mattson,  "LAMMPS  Implementation  of 
Constant  Energy  Dissipative  Particle  Dynamics  (DPD-E) ”,  ARL-TR-XXXX,  U.S.  Army  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD. 
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pair  style  dpde  command 

Syntax: 


pair_style  dpde  kappa_flag  cutoff  seed 


•  kappa_flag  =  0/1  to  turn  off/on  the  energy -dependence  of  kappa  (integer) 

•  cutoff  =  global  cutoff  for  DPD  interactions  (distance  units) 

•  seed  =  random  #  seed  (positive  integer) 


Examples: 


pair_s  tyle 
pair_s  tyle 
pair_coeff 
pair_coeff 


dpde  0  2,5  34387 
dpde  1  2,5  34387 
*  *  3,0  1,0  1,0  1,0 
1  1  3,0  1,0  1,0  1,0 


Description: 

Style  dpde  computes  a  force  field  for  dissipative  particle  dynamics  (DPD)  under  isoenergetic  conditions.  The 
force  on  atom  I  due  to  atom  J  is  given  as  a  sum  of  3  terms: 


f  =  (Fc  +  FD  +  F^r-j 
Fc  =  Aujij 
F°  =  -7<4  (nj  •dTj) 

FR  =  aijuijCij(At)-1/2 


nj  <  rc 


where  Be  is  a  conservative  force,  Pd  is  a  dissipative  force,  and  ft  is  a  random  force.  Rij  is  a  unit  vector  in  the 
direction  Ri  -  Rj,  Vij  is  the  vector  difference  in  velocities  of  the  two  atoms  =  Vi  -  Vj,  zeta  is  a  Gaussian 
random  number  with  zero  mean  and  unit  variance,  and  dt  is  the  timestep  size.  Gamma  is  set  equal  to 
(sigma*  sigma)  I  (2  Kb  Theta),  where  Kb  is  the  Boltzmann  constant  and  Theta  is  the  particle  internal 
temperature. 

For  style  dpde,  the  weighting  factor,  omega_ij ,  varies  between  0  and  1 ,  and  is  chosen  to  have  the  following 
functional  form: 


where  Rc  is  the  cutoff  radius.  Note  that  alternative  definitions  of  the  weighting  function  exist,  but  would  have 
to  be  implemented  with  a  separate  pair  style  command. 

The  kappa_ij  variable  can  be  specified  with  or  without  an  energy  dependence  by  toggling  the  kappa_flag.  In 
the  energy -in  depen  doit  model  (kappa_flag=  0),  kappa_ij  is  explicitly  given  as  a  pair  coefficient.  In  the 
energy -depen  dent  model  (kappa_flag  =  1),  kappa_ij  is  given  by  the  equation: 
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where  kappaO  is  the  pair  coefficient  that  is  specified  in  the  input  file,  kB  is  the  Boltzmann  constant,  and  u_i  is 
the  total  internal  energy  of  particle  I. 

The  differential  internal  conductive  and  mechanical  energies  are  computed  as 


du 


cand 


du 


mech 


•  Vi j)2  -  +  —Wij  -  •  wTi)Ci(Ai) 

2  rtJ  4  mi  mj  2  rl} 


where 


It 

CN 

s 

"IksKij 

•J 

27 ijkB®ij 

0"-1  = 

11  1 
2  0i  +  Oj 

l3 

Zeta_q  is  a  second  Gaussian  random  number  with  zero  mean  and  unit  variance  that  is  used  to  compute  the 
internal  conductive  energy. 


For  style  dpde,  the  pairwise  energy  associated  with  style  dpde  is  only  due  to  the  conservative  force  term  Fc, 
and  is  shifted  to  be  zero  at  the  cutoff  distance  Rc.  The  pairwise  virial  is  calculated  using  only  the  conservative 
term. 


For  style  dpde,  the  following  coefficients  must  be  defined  for  each  pair  of  atom  s  types  via  the  pair  coeff 
command  as  in  the  examples  above,  or  in  the  data  file  or  restart  files  read  by  the  read  data  or  read  restart 
commands: 


•  A  (force  units) 

•  sigma  (force*timeA(l/2)  units) 

•  kappa_ij  (energy *temperatureAime  units)  or  kappaO  (1  /time  units) 

•  cutoff  (distance  units) 

The  last  coefficient  is  optional.  If  not  specified,  the  global  DPD  cutoff  is  used.  Note  that  gamma  is  set  equal  to 
sigma*sigma/(2  Theta),  where  Theta  is  the  average  internal  temperature  of  the  pair. 


Mixing,  shift,  table,  tail  correction,  restart  info: 

The  pair  style  does  not  support  mixing  Thus,  coefficients  for  all  I,J  pairs  must  be  specified  explicitly. 

The  pair  style  does  not  support  thepair  modify  shift  option  for  the  energy  of  the  pair  interaction.  Note  that  as 
discussed  above,  the  energy  due  to  the  conservative  Fc  term  is  already  shifted  to  be  0.0  at  the  cutoff  distance 
Rc. 

pair_style  dpde  command  1 2 
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The  pair  modify  table  option  is  not  relevant  for  this  pair  style. 

The  pair  style  does  not  support  the  pair  modify  tail  option  for  adding  long-range  tail  corrections  to  energy  and 
pressure. 

The  pair  style  writes  its  information  to  binary  restart  files,  so  pair  style  and  pair_coeff  commands  do  not  need 
to  be  specified  in  an  input  script  that  reads  a  restart  file.  Note  that  the  user-specified  random  number  seed  is 
stored  in  the  restart  file,  so  when  a  simulation  is  restarted,  each  processor  will  re-initialize  its  random  number 
generator  the  same  way  it  did  initially.  This  means  the  random  forces  will  be  random,  but  will  not  be  the  same 
as  they  would  have  been  if  the  original  simulation  had  continued  past  the  restart  time. 


Restrictions: 

The  pair  style  dpde  is  only  available  if  LAMMPS  is  built  with  the  USER-DPDE  or 
USER-DPDE-SHARDLOW  package. 

This  pair  style  dpde  requires  the  dpd  atom  style  to  be  used  in  order  to  properly  account  for  the  particle 
internal  energies  and  temperature. 

This  pair  style  dpde  requires  you  to  use  die  communicate  vel  ves  option  so  that  velocites  are  stored  by  ghost 
atoms. 

Related  commands: 

mil'  coeff,  nair  dnd.  fix  dnde 
Default:  none 


(Larentzos)  J.P.  Larentzos,  J.K.  Brennan,  J.D.  Moore,  M.  Lisal  and  W.D.  Mattson,  "Parallel  Implementation 
of  Isothermal  and  Isoenergetic  Dissipative  Particle  Dynamics  Using  Shardlow-Like  Splitting  Algorithms”, 
Submitted  to  Comput.  Phys.  Commun.,  (2013). 


(Larentzos)  I.P.  Larentzos,  J.K.  Brennan,  J.D.  Moore,  and  W.D.  Mattson,  "LAMMPS  Implementation  of 
Constant  Energy  Dissipative  Particle  Dynamics  (DPD-E)”,  ARL-TR-XXXX,  U.S.  Army  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD. 
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pair  style  dpde/conservative  command 

Syntax: 


pair_style  dpde/conservative  kappa_flag  cutoff  seed 

•  kappa_model  =  0/1  to  turn  off/ on  the  energy -dependence  of  kappa  (integer) 

•  cutoff  =  global  cutoff  for  DPD  interactions  (distance  units) 

•  seed  =  random  #  seed  (positive  integer) 


Examples: 


pair_s  tyle 
pair_s  tyle 
pair_coeff 
pair_coeff 


dpde/conservative  0 
dpde/conservative  1 
*  *  3,0  1,0  1 ,0  1,0 
1  1  3,0  1,0  1 ,0  1,0 


2,5  34387 
2,5  34387 


Description: 

Style  dpde/conservative  computes  the  conservative  force  for  dissipative  particle  dynamics  (DPD).  The 
conservative  force  on  atom  I  due  to  atom  J  is  given  by 

Fc  =  Aju'ij  rl}  <  rc 


where  the  weighting  factor  varies  between  0  and  1,  and  is  chosen  to  have  the  following  functional  form: 


a  =  1 


where  Rij  is  a  unit  vector  in  the  direction  Ri  -  Rj,  and  Rc  is  the  cutoff.  Note  that  alternative  definitions  of  the 
weighting  function  exist,  but  would  have  to  be  implemented  with  a  separate  pair  style  command. 


The  kappajij  variable  can  be  specified  with  or  without  an  energy  dependence  by  toggling  the  kappa_flag.  In 
the  energy -in  depen  doit  model  (kappa_flag=  0),  kappa_ij  is  explicitly  given  as  a  pair  coefficient.  In  the 
energy -depen  dent  model  (kappa_flag  =  1),  kappajj  is  given  by  the  equation: 

Ui  +  Uj 

2 


- 

kB 


where  kappaO  is  the  pair  coefficient  that  is  specified  in  the  input  file,  kB  is  the  Boltzmann  constant,  and  u_i  is 
the  total  internal  energy  of  particle  I. 

This  pair  style  differs  from  the  other  dpd  styles  in  that  the  dissipative  and  random  forces  are  not  computed 
within  the  pair  style.  This  style  is  combined  with  the  fix  dpde/shardlow ,  which  will  perform  the  stochastic 
integration  of  the  dissipative  and  random  forces  through  theShardlow  splitting  algorithm  approach. 


For  style  dpde/conservative,  the  pairwise  energy  associated  with  style  dpde/conservative  is  only  due  to  the 
conservative  force  term  Pc,  and  is  shifted  to  be  zero  at  the  cutoff  distance  Rc.  The  pairwise  virial  is  calculated 
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using  only  the  conservative  term. 

For  style  dpde/consen’ative,  the  following  coefficients  must  be  defined  for  each  pair  of  atoms  types  via  the 
pair  coeff  command  as  in  the  examples  above,  or  in  the  data  file  or  restart  files  read  by  the  read  data  or 
read  restart  commands: 

•  A  (force  units) 

•  sigma  (force*timeA(l/2)  units) 

•  kappa  ij  (energy* Temperature/time  units)  or  kappaO  (1/time  units) 

•  cutoff  (distance  imits) 

The  last  coefficient  is  optional.  If  not  specified,  the  global  DPD  cutoff  is  used.  Note  that  gamma  is  set  equal  to 
sigma*sigma/(2  Theta),  where  Theta  is  the  average  internal  temperature  for  the  pair. 


Mixing,  shift,  table,  tail  correction,  restart  info: 

The  pair  style  does  not  support  mixing.  Thus,  coefficients  for  all  IJ  pairs  must  be  specified  explicitly. 

The  pair  style  does  not  support  the  nair  modify  shift  option  for  the  energy  of  the  pair  interaction.  Note  that  as 
discussed  above,  the  energy  due  to  the  conservative  Fc  term  is  already  shifted  to  be  0.0  at  the  cutoff  distance 
Rc. 

The  pair  modify  table  option  is  not  relevant  for  these  pair  styles. 

The  pair  style  does  not  support  the  nair  modify  tail  option  for  adding  long-range  tail  corrections  to  energy  and 
pressure. 

The  pair  style  writes  its  information  to  binary  restart  files,  so  pair_style  and  pair_coeff  commands  do  not  need 
to  be  specified  in  an  input  script  that  reads  a  restart  file.  Note  that  the  user-specified  random  number  seed  is 
stored  in  the  restart  file,  so  when  a  simulation  is  restarted,  each  processor  will  re-initialize  its  random  number 
generator  the  same  way  it  did  initially.  This  means  the  random  forces  will  be  random,  but  will  not  be  the  same 
as  they  would  have  been  if  the  original  simulation  had  continued  past  the  restart  time. 


Restrictions: 

The  pair  style  dpde/conservative  is  only  available  if  LAMMPS  is  built  with  the  USER-DPDE-SHARDLOW 
package. 

The  pair  style  dpde/conservative  requires  the  dpd  atom  style  to  be  used  in  order  to  properly  account  for  the 
particle  internal  energies  and  temperature. 

The  pair  style  dpde/conservative  requires  you  to  use  the  communicate  vej  yes  option  so  that  velocites  are 
stored  by  ghost  atoms. 

The  pair  style  dpde/consetvative  will  not  restart  exactly  when  using  the  read  restart  command,  though  they 
should  provide  statistically  similar  results.  This  is  because  the  forces  they  compute  depend  on  atom  velocities. 
See  the  read  restart  command  for  more  details. 

Related  commands: 
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pair  coeff.  pair  dpd.  pair  dpde.  fix  dpde/shardlow 


Default:  none 


(Larentzos)  J.P.  Larentzos,  J.K.  Brennan,  J.D.  Moore,  M.  Lisal  and  W.D.  Mattson,  "Parallel  Implementation 
of  Isothermal  and  Isoenergetic  Dissipative  Particle  Dynamics  Using  Shardlow-Like  Splitting  Algorithms”, 
Submitted  to  Comput.  Phys,  Commun.,  (2013). 


(Larentzos)  J.P.  Larentzos,  J.K.  Brennan,  J.D.  Moore,  and  W.D.  Mattson,  "LAMMPS  Implementation  of 
Constant  Energy  Dissipative  Particle  Dynamics  (DPD-E)",  ARL-TR-XXXX,  U.S.  Army  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD. 


pair_style  dpde/conservative  command 


16 


19 


1  DEFENSE  TECHNICAL 
(PDF)  INFORMATION  CTR 

DTIC  OCA 

2  DIRECTOR 

(PDF)  US  ARMY  RESEARCH  LAB 
RDRL  CIO  LL 

IMAL  HRA  MAIL  &  RECORDS  MGMT 

1  GOVT  PRINTG  OFC 
(PDF)  AMALHOTRA 

1  DIR  USARL 
(PDF)  RDRL  WML  B 
J  LARENTZOS 


20 


